Temperature-induced crossovers in the static roughness of a one-dimensional interface 
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At finite temperature and in presence of disorder, a one-dimensional elastic interface displays 
different scaling regimes at small and large lengthscales. Using a replica approach and a Gaus- 
sian Variational Method (GVM), we explore the consequences of a finite interface width £ on the 
small-lengthscale fluctuations. We compute analytically the static roughness B(r) of the interface 
as a function of the distance r between two points on the interface. We focus on the case of short- 
range elasticity and random-bond disorder. We show that for a finite width £ two temperature 
regimes exist. At low temperature, the expected thermal and random-manifold regimes, respec- 
tively for small and large scales, connect via an intermediate 'modified' Larkin regime, that we 
determine. This regime ends at a temperature-independent characteristic 'Larkin' length. Above 
a certain 'critical' temperature that we identify, this intermediate regime disappears. The thermal 
and random-manifold regimes connect at a single crossover lengthscale, that we compute. This is 
also the expected behavior for zero width. Using a directed polymer description, we also study via a 
second GVM procedure and generic scaling arguments, a modified toy model that provides further 
insights on this crossover. We discuss the relevance of the two GVM procedures for the roughness at 
large lengthscale in those regimes. In particular we analyze the scaling of the temperature-dependent 
prefactor in the roughness B{r) ~ T 2]p r 2( and its corresponding exponent p- We briefly discuss the 
consequences of those results for the quasistatic creep law of a driven interface, in connection with 
previous experimental and numerical studies. 



I. INTRODUCTION 

Almost everyone has already unwittingly spilt coffee 
on her/his work table and observed the inexorable pro- 
gression of the liquid into her/his favorite article. How- 
ever inconvenient may be the consequences of such a sim- 
ple home experiment, it shares in fact a lot of common 
physical features with a variety of systems and phenom- 
ena ranging from ferromagnetic 1 ' 2 and ferroelectric 3, do- 
main walls, to growth surfaces 5 , contact line in wetting 
experiments 6,7 , or crack propagation in paper 8 . All those 
systems display different coexisting phases separated by 
an interface, whose shape and dynamics are determined 
from two competing tendencies: the elastic cost of the 
interface which tends to flatten it while disorder in the 
environment induces deformations adapting the interface 
to the local energetic valley and hills. 

To describe these phenomena, a successful theoretical 
approach is that of disordered elastic systems 9,10 (DES), 
in which the bulk details of the interface are summa- 
rized in its mere position, seen as a fluctuating man- 
ifold whose energy is the sum of elastic and disorder 
contributions. Such a description also encompasses pe- 
riodic systems such as charge density waves 11,12 or vor- 
tex lattices 13 occurring in type II superconductors. This 
approach accounts for a complex free-energy landscape 
which exhibits metastability and explains the dynamical 
glassy properties - such as hysteresis, creep or ageing - 
observed in experimental realizations. Such systems dis- 
play similar features independently of the scale at which 
they are observed: in other words, they present ranges 
of scale on which they are statistically scale-invariant 14 . 



A simple way to characterize this property is to study 
the roughness of the interface (denoted B(r)), defined as 
the variance of the relative displacements of the manifold 
at a given lengthscale r. Scale invariance translates into 
having the roughness behaving as a power law B(r) ~ r 2l > 
characterized by a roughness exponent £. 

In addition to its direct experimental relevance e.g. 
in ferromagnetic domain walls, the ID interface shares 
many universal features with other distinct physical 
problems, such as the directed polymer (DP) in random 
media (for reviews, see Refs. 15 and 16), the noisy Burg- 
ers' equation in hydrodynamics 17 ' 18 , or bosons with at- 
tractive hard-core interaction in one dimension 19 . These 
systems are all falling in the so-called Kardar-Parisi- 
Zhang (KPZ) class 20 ' 21 . This problem has now been 
studied for many decades, the attention being focused on 
the large-scale - or so-called 'random-manifold' (RM) - 
properties of the roughness: many approaches have been 
used to establish the non-trivial exponent Crm = 2/3 for 
B(r) ~ r 2 ^ in the large r limit, ranging from dynamical 
renormalization group 17 ' 22 , to hidden symmetries 23 and 
Bethe Ansatz computations 19 . 

In spite of the versatility of the previous studies, there 
has been a recent uprise of interest in different direc- 
tions: from a mathematical point of view it is only very 
recently 24 that the 2/3 RM exponent has been proven, 
the method itself making the link with another class of 
systems - transport models known as asymmetric ex- 
clusion processes. Besides, in a recent series of works 
in the mathematics 25,26 and physics 27 ' 28 communities, it 
has been shown that the free-energy distribution for the 
polymer endpoint is obtained from the Tracy- Widom dis- 
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tribution, when the random potential is delta-correlated. 
Similarly, another interesting and related question which 
seems to have been neglected up to recently is the role of 
temperature in the large scale behaviour of the rough- 
ness: although physically one could a priori expect 
the roughness to be temperature-independent at large 
lengthscales, since it is dominated by the disorder, it hap- 
pens that scaling arguments show this shall not be the 
case at small enough disorder correlation length £, and 
that there is a subtle interplay between low temperature 
and small £ limits 28 ' 29 . 

Despite these recent progress several questions remain. 
In particular, in these works, the width of the interface 
or the disorder correlation length £ is assumed to be zero. 
It is thus important to ascertain how keeping £ finite - as 
always the case in physical realizations - influences the 
physics at stake. The role of £ in the large scale limit is 
exemplified specially in the context of functional renor- 
malization group (FRG) 30 ' 31 , where at large scale the 
roughness arises from the zero-temperature fixed point 
of the renormalization flow, hence displaying no temper- 
ature dependence - in contrast to pure scaling predic- 
tions (in other words, having £ finite seems to influence 
temperature scaling exponents at large distance 29 ). 

In this paper, we address the issue of the role of £ 
concerning the roughness scaling and the temperature 
regimes. We focus on the roughness of a one-dimensional 
static interface subjected to 'random-bond' disorder, at 
finite temperature. It is known that at short lengthscale 
disorder plays no role and the interface is in a thermal 
regime (Cth = |), wnne at large lengthscale the interface 
is in the RM regime, characterized by a roughness expo- 
nent Crm = § (see Refs. 17, 22, and 23). We examine in 
detail the way the roughness crosses over from thermal to 
RM behaviour, and in particular whether there is one or 
several crossover lengthscales, possibly allowing for a non 
trivial 'intermediate regime'. We tackle these problems 
by using a Gaussian variational method (GVM) 32 ~ 35 . Al- 
though this method is only approximate, it allows rather 
complete calculations for the roughness B(r). The al- 
ternative methods are not so well suited to analyse this 
issue for a one dimensional interface. Even though nu- 
merical methods such as Monte-Carlo 36 , Langevin 37 or 
transfer matrix 38 methods are efficient to deal with the 
interface, one would need to implement them for very 
large system sizes, in order to tackle the issue of crossover 
lengthscales. As for the functional renormalization group 
(FRG) approach 30,31 , it is based on an expansion in 
e = 4 — d for a manifold of internal dimension d and 
may not be suited for our case of interest (d = 1). 

The plan of the paper is as follows. In Sec. II we in- 
troduce the DES description of a ID interface. We show 
how the replica trick allows to make the average over 
disorder. In Sec. Ill we detail the variational procedure, 
originally introduced by Mezard and Parisi 33 , that we 
use in the paper to obtain an approximation of the vari- 
ous physical properties of the interface, and in particular 
the correlation function B(r). The corresponding rough- 



ness is computed in Sec. IV, along with the crossover 
lengthscales separating its different power-law regimes, 
and their temperature dependence. In Sec. V we present 
a second GVM procedure, this time in direct-space rep- 
resentation, on a 'toy model' which has been argued to 
be an effective model for the study of a directed poly- 
mer's end-point fluctuations, and on which the model of 
a ID interface can carefully be mapped. In Sec. VI we 
examine how generic scaling arguments shed light on the 
interface and the toy- model results. The two GVM re- 
sults are compared and discussed in Sec. VII, and we 
finally conclude in Sec. VIII. We present for reference 
standard Flory arguments in Appendix A, and we recall 
in Appendix B some useful properties of hierarchical ma- 
trices. 



II. MODEL OF A ID INTERFACE 
A. Model 

In the DES framework, an interface can be described 
as an elastic manifold of dimension d with m transverse 
components, submitted to the random potential of a 
physical space of dimension D = d + m. The coordi- 
nates in this space can thus be split between the internal 
and transverse coordinates of the interface, respectively 
(z, x) e R d x W n = R D . In the case d = m = 1 (simply 
denoted '1+1' for DPs), a one-dimensional interface is 
thus described as an elastic line living in a bidimensional 
plane, with a disordered energy landscape that will be 
defined below. Restricting ourselves to the case without 
bubbles nor overhangs, each configuration of the interface 
can be indexed by a univalued displacement field u(z) 
which parametrizes the position (z,u(z)) of the interface 
along its internal direction z. This is schematically indi- 
cated in Fig. 1(a). 

In disordered systems there is no strict spatial trans- 
lational invariance, but it can be recovered statistically 
once the disorder is averaged out. Thereafter, we will 
thus mostly work in Fourier-transform representation 
and denote the Fourier modes along z and x respec- 
tively q and A. However there are in fact both infra-red 
and ultra-violet cutoffs in those Fourier modes, since a 
physical realization of an interface lives in a finite bidi- 
mensional plane of typical size L and is supported by a 
sublattice (e.g. a crystal in a solid) whose spacing 1/A 
defines the smallest physical lengthscale in the system. 
The cutoffs 1/L < q, A<A will be conveniently reintro- 
duced whenever needed to cure non-physical divergences 
in Fourier-transform integrals. Thereafter we use the no- 
tation 3 d q = jj^fa, as well as 8^ d \q) = (2ir) d 5^(q), and 
similarly for other Fourier modes. 

More importantly, in order to implement a finite inter- 
face width £ into this model, we introduce the density of 
the interface p u [z,x) which is positive and normalized 
at fixed z by J R dx ■ p u (z, x) = 1, decreases significantly 
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We can now construct the DES Hamiltonian of the 
interface, which is the sum of the energetic cost of its 
distortions H e i [u] and the contribution of the disordered 
energy landscape Hdis [u,V]. We assume that the elas- 
tic limit |V z it(z)| <^C 1 is realized and that the elastic- 
ity is short range, so that the energy per Fourier mode 
u q = j R dz ■ u{z)e lqz is cq 2 , where c is the elastic con- 
stant. The mode u q= o corresponds to the mean position 
of the interface, and introduces an additive constant in 
H c i which will disappear in the Boltzmann weight and 
moreover does not contribute to the roughness. It can 
thus be put equal to zero directly in the elastic Hamil- 
tonian by the ad hoc redefinition of u. We assume a 
random-bond disorder, i.e. that the interface couples lo- 
cally to the random potential. Thus the full Hamiltonian 
% is given by: 

H [u, V] = U ei [u] + H dis [u, V] (3) 
%a M = | J dz ■ (V z w(z)) 2 = * J dq-cq 2 U- q u, 



Hdis [u, V] — / dz dx ■ p u (z, x) V(z, x) 
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(4) 
(5) 



Note the alternative formulation of the Hdis with an 
effective random potential V coupled to a zero-width in- 
terface: 



u,V = dz ■ V (z, u{z)) 



Hdis 



and its associated disorder distribution 



(6) 



(b) 



FIG. 1. (a) Displacement from a given flat configuration de- 
scribed by the z axis, (b) Gaussian density p u (z,x) centered 
on the position (z,u(z)) of the interface. 



for | at — u(z)\ > £ and tends to a Dirac ^-function when 
£ — > 0. As schematized in Fig. 1(b) we choose the sim- 
plest smooth density which is a Gaussian centered on 
(z,u(z)), of uniform standard deviation £, whose Fourier 
representation is thus: 



Pu (z,x)= / d\ ■ e ^-u(z)) e -x*e,2 



(1) 



The interface is subject to a random potential V{z, x) 
whose distribution is Gaussian and uncorrelated in space. 
This corresponds to the limit of the collective pinning by 
many weak impurities. Denoting the statistical average 
over disorder by an overline, since it is Gaussian it is fully 
defined by: 



V(z,x) = 0, V(z,x)V(z',x') = D-5{z-z')5{x-x') (2) 
where D is the strength of disorder. 



V(z, x) = 0, V(z, x) Viz', x') = 6(z - z') R 6 (x - x') (7) 

where R^(u) is up to an additive constant the usual cor- 
relator of the disorder which is the key function renor- 
malized in FRG procedures on DES 31 ' 34 . It corresponds 
to the overlap between two densities p u in our formu- 
lation and encodes the statistical translational invari- 
ance of the random potential. Note that R^(u z ,u' z ) = 
D J R dx ■ p u (z,x)p u '(z,x) inherits both the translational 
invariance and the symmetry of the density. The Gaus- 
sian density (1) translates into the following Gaussian 
correlator: 



D 



d\ ■ 



(8) 



So the parameter £ can be seen either as the width 
of the interface or as the disorder correlation length (or 
a convolution of both). In previous GVM computations 
on DES 33 ' 34 the correlator R{u) was assumed to exhibit 
an asymptotic power-law behaviour whose exponent de- 
pended on the university class of the disorder. Here we 
focus on the role of a finite £ rather than on this asymp- 
totic behaviour, following closely a similar approach of 
periodic DES 32 . 
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B. Statistical averages 

The static properties of an interface are accessible by 
averaging over its thermal fluctuations and the stochas- 
tic variable V which is associated to each configuration 
of quenched disorder, so two statistical averages have to 
be successively performed for a given observable O. The 
first one is the thermal average at fixed disorder (0) v us- 
ing the Boltzmann weight e~ m[uy] / Z v , where /3 = 1/T 
is the inverse of the temperature (taking the Boltzmann 
constant ks — 1) and Zy the canonical partition func- 
tion: 

Z V = Jvu- e -m*y\ (9) 

(0) = J- [vu-0[u]- e-WlvM (io) 
Zv J 

where the functional integral J T>u sums over all possible 
configurations u of the interface. 

The second one is the average over disorder O which 
has already been introduced in (2), assuming that the 
system is large enough to be disorder self-averaging (the 
equivalent of ergodicity for thermal averages) 35 . To re- 
cover a translational invariance and be able to work in 
Fourier space, one would like to average technically first 
over disorder. This can in fact be done using the well 
known replica trick 35 . Indeed, introducing n replicas of 
the partition function at fixed disorder Zv, n being an 
arbitrary integer, we can replace 1/Zy in the thermal av- 
erage (10) by lim n _>o Zy n ~ x under the (strong) assump- 
tion that at the end of our computations the analytical 
continuation n is well-defined and physically mean- 
ingful. This gives: 

(0) v = lim / Vu 1 {---)Vu n -0[u 1 ] . e -^E: =1 «K,v] 

n-+QJ 

(11) 

Using the linearity in V of Hdis (5) and the Gaussian dis- 
tribution of disorder (2), we can then explicitly average 
over disorder and define an effective replicated Hamilto- 
nian % which couples all the replicas (u = (ui, ...,«„)): 



(O) = lim / V Ul {- ■ ■ )Vu n ■ O [u{\ ■ e -^ = i«K^] 

TWO J 

= lim / Vui{- ■ ■ )Vu n ■ O [ui] ■ e - 0il W 



n-yO 



(12) 



where 



n {3 f U 

M-M = ^ ^el [u a ] - 2 dz ■ ^ R ( (u a (z) - u b ( 

o=l ^ R a, 6=1 



(13) 

We have thus reformulated the problem of one interface 
in a random potential V into a system of n coupled in- 
terfaces without disorder, in the limit n — > 0. 



Finally, using (4) and (8) , we explicit the effective repli- 
cated Hamiltonian which is exact up to this point, diag- 
onal in Fourier space due to translational invariance and 
depends only on the parameters {£, c, D, T}: 



H [u] = H c \ [u] + Hdis [u\ 

If " 

^oi [u] = - / dq- cq 2 V u a (-q)u a (q) 

2 Jr a=l 
U dia ffl=-iR f d\-e~^ 2 fdz-f" e^^W— 



(*)) 



,6=1 



(14) 



Roughness and displacement correlation 
function 



In order to study the static fluctuations of the position 
of the interface, we compute the variance of the relative 
displacements of two points of the interface: B(zi, Z2) = 



(u(zi) — u{z2)) )■ The disorder average leads back to 

translational invariance B{zx,z%) = B{z\ — ^2,0). One 
can thus define the roughness as a function of the length- 
scale r, formally the two-point correlation function of our 
system: 



(15) 



which is the Fourier transform of the structure factor 
S(q): 



S{q) = / dq - (u-qU q ) 



B{r) = J dq ■ 2 (1 - cos (qr)) S(q) 



(16) 



If the DES displays a scale invariance in its fluctua- 
tions, the roughness is expected to follow a power-law 
behaviour B{r) ~ r 2l », with a corresponding roughness 
exponent Q. To describe the possible interplay between 
different terms in B{r) we generalize the definition of £: 



1 d\og B(r) 
^ [r> ~2 <91ogr 



(17) 



whose different values characterize the different regimes 
of fluctuations depending on the lengthscale considered. 
Along with the r-independent prefactor of B(r) it probes 
the physics at different lengthscales and a roughness 
regime is actually defined by a constant value of (. 

If the exact replicated Hamiltonian could be put into 
a quadratic replicated form, diagonal in Fourier space, 
such as: 



If " 



(18) 



a, 6=1 
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then the corresponding structure factor would be: 

S (q) = /T 1 • Urn G aa (q) 
n—>0 



(19) 



In the absence of disorder (H = Hd) the replicas are 
uncoupled and G~ b (q) = cq 2 ■ 8 ab , so the structure fac- 
tor is proportional to the 'thermal propagator' -K which 

cq 

leads to the purely thermal roughness of a ID inter- 
face B t h(r) = —■ The exact Hamiltonian (14) is di- 
agonal in Fourier space, but it cannot be put into a 
quadratic form because of the coupling of all replicas 
E",6=i e lA( "' (z) "" i(z)) - This forbids the direct use of (19) 
for the computation of S(q). To be able to do so, we 
thus approximate % in the Boltzmann weight (12) by a 
quadratic replicated Hamiltonian Ho optimized by GVM, 
as explained in the next section. 



III. GVM AND FULL-RSB ANSATZ 

The GVM has already been applied to DES, periodic 
systems 32 as well as manifolds 33,34 , to study the temper- 
ature dependence of observables at thermodynamic equi- 
librium. Here we extend these computations specifically 
to the case of a one-dimensional interface of finite width 
£, in order to explore its small-lengthscales behaviour. 

We follow in this derivation the main steps outlined in 
Ref. 32 for periodic systems. One important difference 
comes from the fact that in our case the variable A in (14) 
is continuous, while it takes discrete values for periodic 
systems. As we will see this has drastic consequences for 
the physical properties of the system as well as for the 
calculation itself. 



A. GVM in Fourier representation 

The variational method consists in replacing, in the 
Boltzmann weight of statistical averages, the exact 
Hamiltonian % or more generally the exact action of a 
system by a trial Hamiltonian Hq with variational param- 
eters. The criterion chosen to optimize this approxima- 
tion is given by the Gibbs-Bogoliubov inequality, which 
states that the free energy J 7 of a system is minimum at 
equilibrium, i.e. when the probability measure of the sys- 
tem is precisely described by its exact Boltzmann weight: 



(20) 



where F va , T is the variational free energy associated to the 
trial Hamiltonian T-Lq and which has to be minimized, T 
and J-"o the free energies corresponding respectively to 
H and %q (they are defined with respect to their cor- 
responding partition function Z = e~ ,3 ~ F ), and (O) the 
statistical average defined over Ho- 

For the replicated Hamiltonian (14), the trial Hamil- 
tonian 'Hq is chosen quadratic of the generic form (18), 



parametrized 33 by a n x n matrix G~ b (q). Its inverse 
matrix gives thus directly access to the correlation func- 
tions: 



(u a (-q) Ub(q)) = 13 1 • G ab (q) ■ S (q - q) 



(21) 



and in particular to the structure factor S(q) and the 
roughness B{r) itself. The minimization of J- vax with 
respect to the variational parameters G a b(q) (the Green 
function of "Ho) gives a saddle point equation for the opti- 
mal matrix G~ b (q). Besides, the replica trick constrains 
the structure of G~ b , which must be a hierarchical ma- 
trix. In Appendix B we recall some useful properties of 
such matrices, including their inversion formulas in the 
limit n — > which will be used extensively thereafter. 

The extremalization condition dT var /dG a b{q) = can 
be reformulated as: 



d 



&ab 



8 0G oh ( q ) {: ' 



(22) 



Note that a ab is independent of the Fourier mode q since 
Hdis is itself purely local in z. 

We point out that although the GVM approach is 
only approximate in our context, it becomes exact 34 in 
the limit of an infinite number of transverse components 
772 — )■ oo. The extremalization equations then appear as 
genuine saddle-point equations with 1/m playing the role 
of a small parameter. By extension we extensively use 
thereafter this improper denomination to refer to (22). 
Furthermore, one could for completeness check the sta- 
bility of the GVM solution by considering its associated 
Hessian matrix 33 . This problem can be quite compli- 
cated, so we simply check the physical consistency of our 
results in what follows. 



B. Saddle point equation in the full-RSB 
formulation 

To compute explicitly a a b (22) we start from (14) and 
performing the Gaussian statistical average 



iX(u a (z)-u b (z)\ _ e -^-{(uJz)-u b (z)) 2 ) Q 



(23) 



we have: 



W, ) - PD 



[ dz I d\-J2 e -A a tt a +$««.W-«>W>">o> 
Jr Jr b=1 



(24) 

where the variance of the relative displacement between 
two replicas at fixed z is obtained by applying (21) 
on its Fourier-transform representation. This gives the 
translational- invariant quantity: 



(K(z) - u b {z)f) =T dq (G aa (q) + G bb {q) - 2G ab {q)) 
\ '0 Jr 

(25) 
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Performing d/dG a b{q) and denoting G aa = G we even- 
tually obtain a a b, separating off-diagonal a b terms 

a ^ b = V f dX • A 2 • e - x2 ^ +T i< d ^- G ^)) (26) 



from the diagonal a = b: 



(27) 



Note that {Tidisj is extensive in the interface size 

(see (24)) whereas a a b is intensive as expected since J dz 
has disappeared. 

We can compute the connected part of the hierarchical 
matrix G^ 1 (q): 

G-\q) = £ G£,(g) = eg 2 - a - £ ov^a = cq 2 (28) 



Using (B3) we know that: 
G c (q) 



cq- 



(29) 



In this GVM framework the connected part of the hierar- 
chical matrices G~ 1 (q) and G(q) have thus a straightfor- 
ward physical meaning: in the Hamiltonian it represents 
the elastic energy per Fourier mode G~ 1 {q) = cq 2 and in 
the Green function it gives back the thermal propagator 
Gc(q) — ^2 • On one hand the irruption of disorder popu- 

lates by construction the off-diagonal elements of G^ 1 (q) 
with the ^-independent coupling terms —a a ^b under the 
constraint (28). On the other hand the invariance of 
G c {q) is a consequence of the statistical tilt symmetry of 
the DES description 31 ' 32 ^ 39 - 41 . 

To determine the actual propagators we consider the 
generic case of a full replica-symmetry breaking (RSB) 
Ansatz, the off-diagonal term being parametrized by u € 
[0,1]: 

G-\q) = (G-\q) - or, -a(u)) G(q) = (d(q), G(q, u) 

(30) 

where the definition of the connected part G c 1 is consis- 
tent with the continuous version (B8) for the definition 
of a (27) . The saddle point equation (26) becomes: 



a(u) 



D 
T 

D 
4ttT 



cf A • A ■ ■ 



-A 2 (« 2 +T / <tq(G(q)-G(q,u))) 



(31) 



e 



T 



dq(G(q)-G(q,u) 



-3/2 



(32) 



and the relation between (^G(q) — G(q,u)j and cr(u) can 

be made explicit using the inversion formulas (Bll) and 
(B12) and the definition (B9) of the self-energy [a] (u) 
which appears in them. a(u) must thus satisfy this self- 
consistent saddle-point equation, but we will work alter- 
natively with er(0) and [a] (u). 



C. Determination of tr(lt) and [a] (u) 

We assume that a(u) is continuous by sector, and aim 
to point out its power-law behaviour using the following 
implication of the definition (B9): 



[a] (u) = u ■ cr'(u) 



(33) 



We first apply d u on (32), then identify a power of a(u) 
and use (B12) and (B15) to make explicit the following 
term: 



T J dq [G(q)-G(q,u 



T dq-d u 

JR 
T f dq 



. u (cq 2 + [a}(v)) 2 
a'(u) T 



(eg 2 + [a] (u)) 2 4^ 



We thus obtain the new equation: 



a'(u) ■ [a] (u)- 3 / 2 
(34) 



37T 1 / 3 



a'(u)=a'(u).(Ta(u)f 3 W (u)^ 

(35) 

We have either u'{u) = which corresponds to plateaux 
in a{u) and possibly to a full replica-symmetric (RS) so- 
lution, or if <t'(u) 7^ we have a strictly monotonous 
segment that should satisfy the new equation: 



(36) 



Note that by definition of the self-energy (B9) we have 
[a] (0) = 0, so the relation (36) implies in particular that 
ct(0) = 0, property that can be checked a posteriori by 
plugging our full-RSB solution directly into the initial 
saddle point equation (32). 

Differentiating again (36) and using (33) to reintroduce 
a u dependence, we finally obtain for the monotonous 
segments of the self-energy [a] (it): 



[a] (u) 



,14 



c 3 ^ 4 • ( U /T) 10 ^ (C , D) • ( U /T) iu (37) 
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b W TT 2 



It is important to note that the power-law form (u/T) , 
which will actually condition the asymptotic tempera- 
ture dependence of the roughness, is totally constrained 
by the GVM procedure. Indeed, we can track down the 
temperature factors in the previous procedure: the ex- 
tremalization condition of Jvar introduces a first j3 in (22) 
via dJ^o/dGabiq) — (3~ 1 G~ b (q). The derivative involv- 
ing ^"Hdis^ introduces no new T-dependence because 
its overall f3 factor is canceled through the derivative 
d/dG a b(q) of its exponential argument. The only T fac- 
tor which remains in the relation (36) between a(u) and 
[a] (u) eventually leads to the (u/T) u dependence once 
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the u dependence is made explicit. We will go back to 
that property in Sec. IV E, and compare its implications 
for the asymptotic T-dependence of the roughness for 
the ID interface versus the prediction of another model 
in Sec. VII. 

We have obtained so far that for the values of u for 
which cr(u) is not constant, it must obey 

10 A (c ^ D) 



cr(u) 



(u/TY 



cte 



(38) 



9 T 

where this last constant has to be determined in rela- 
tion to the possible plateaux and cutoffs imposed on the 
Ansatz er(u), using the saddle point equation (32) to fix 
them. 

In addition to the above power-law behaviour we must 
consider the possibility of a plateau in a(u). One such 
trivial example would be a totally replica-symmetric 
(RS) solution for which cr(u) is independent of u. In 
that case, in the saddle-point equation (31) the term 

^2 as explained in (B5) and 



G(q)~G(q,u) = G c {q) = 
thus: 

D f 

ct R s(u) = ctrs(0) = — / 
1 Jr 



d\-X 2 -e- x ^ 2 e- x2 S^ G ^ =0 



(39) 

The RS solution completely eliminates all effects of dis- 
order, it is thus clearly unphysical, as was the case in the 
higher dimensional cases 33 and the periodic ones 32 . 

To take into account the possible presence of plateaux 
we thus look for a full-RSB solution cr(u) with a single 
cutoff v c £ (0,1) of the form: 

10 A {C}D) 



a (u < v r ) 

9 T 

a (u > v c ) = a (v c ) 



{u/Tf+a{Q) 



(40) 



(41) 



and for the self-energy: 

M (u<v c )=A (c>D y(u/T) 10 
[a] (u > v c ) = [a] (v c ) 

Up to this point the width £ does not appear in these 
expressions. It is in fact encoded into the single full- 
RSB cutoff u c (£). The equation for this cutoff is ob- 
tained by checking the consistency of the solution (41) 
with the initial saddle point equation (32). Using the in- 
version formula (B12) adapted to the presence of a cutoff 
a' (u > v c ) = 0, then integrating over the q modes with 
(B14) and (B15) and inserting our full-RSB solution (41) 
we have: 



T 



I dq ■ (G(q) 
Js. v 



T 



[<r] (v c ) 



-1/2 



2,/c 



G(q,u) 



T 



dv ■ cr'(v) 



[a] (v) 



-3/2 



~ * • F 

2-3 7 (cDY 



6 VT 



(42) 



Substituting this expression into (32), the full-RSB 
Ansatz u (u < v c ) (40) is self-consistent if the in- 
dependence cancels the v c terms: 



5° 1 

:1T- 



2-3 7 ( cD y VT 



(« c - 5/6) = (43) 



and also if <r(0) = 0. This last condition is enforced by 
(36) and can be checked a posteriori by combining (32), 
(42) and (43): 



o-(O) = — / d\ ■ X A ■ lim e" 



(44) 



The single full-RSB cutoff is thus given by the polynomial 
equation: 



vt = A (5/6 -v c ), A 



5 5 T 6 
:7T- 



2-3 7 {icD f 



(45) 



whose solution is plotted in Fig. 2(b). The factor 5/6 
can be shown to be closely related to the Flory exponent 
(f in d — 1 (cf. Appendix A); the presence of a cutoff 
v c < 5/6 = (2CF|rf=m=i)~ 1 in the full-RSB Ansatz [a] (u) 
is in fact necessary for [a] (it) to be a solution of the GVM 
saddle point equation. 

All the above results are summarized in Fig. 2. 



D. Low versus high temperature limits of [a] (u) 

Before computing the roughness B(r), which we will 
do in the next section, we have to explicit the low and 
the high temperature limits of this solution, since we also 
aim to probe the T-dependence of the roughness. 

An explicit analytical expression for v c € (0,5/6) can 
be obtained for the two opposite limits of A, which at 
fixed £, D > correspond respectively to T — > and 
T — > oo and yield the following asymptotic behaviour: 



A^o 5 /4 _ 

v c « - 

A— >oo , 

v c w 5/6 + 0" 



1/3 



T 



(Cc£>) 1 /3 



(46) 
(47) 



The crossover between the two regimes of high and 
low temperature is in fact conditioned by the value of 
the dimensionless parameter A in the equation for v c (£) 
(45); an arbitrary definition of a 'critical' temperature is 
naturally given by the crossover between the two opposite 
limits A — > oo and A — > 0, which happens at: 



A = l 



Tr 



{&Dy/ 3 



2-3 7 

5 5 7T 



1/6 



0.87 



(48) 



and this last constant is of order 1. Note that this crit- 
ical temperature depends explicitly on the width £, in 
such a way that the limits T — > and £ —> can- 
not be exchanged with impunity: imposing £ = from 




thermal fluctuations: 



(a) 



5/6 -- 








1 



(b) 



FIG. 2. GVM solution for the ID interface, (a) Self-energy 
[a] (u) given by (41); using [a] (it) = ua'(u) and er(0) = 
the expression (40) for a(u) can be recovered, (b) Full-RSB 
cutoff v c as a function of A ~ {^cD)' A ' obtained by solving 
the polynomial equation (45). It starts linearly at A — > 
and saturates to 5/6 at yl — > oo, yielding in particular the 
low-temperature dependence given by (46). 



the beginning is equivalent to considering exclusively 
the 'high' temperature regime, and the somehow non- 
physical regime at simultaneously zero-temperature and 
zero-width has to be carefully handled. 

The first limit (46) implies that when the thermal fluc- 
tuations are suppressed, [a] (u) tends to a non-zero RS 
solution, since v r and 



M («c) 



A->0 1 



4 VVt? 



2/3 



r 10/3 c -l/3^2/3 (49) 



Increasing the strength of disorder D is compatible with 
this limit, and indeed the self-energy is intuitively ex- 
pected to increase with D. 

On the contrary the decrease of D or £ leads to the 
second limit (47), in which the ^-dependence has been 
washed out from the GVM solution by the relatively large 



M («c) 



3 4 c 3 D 4 
2*% T 10 



(50) 



IV. GVM ROUGHNESS AND CROSSOVER 
LENGTHSCALES OF THE ID INTERFACE 

Using the definition of the structure factor (16) in rela- 
tion with the Green function of a quadratic Hamiltonian 
(19), we can now compute the corresponding roughness 
as a function of the lengthscale r along the internal co- 
ordinate z of the interface: 



B(r) = T 



2(1 - cos(qr)) ■ lim G(q) 

n— s-0 



(51) 



The inversion formula for lim n _>o G(q) (BIO) contains 
two contributions (since er(0) = 0) which yield respec- 
tively a purely thermal and a disorder-induced roughness: 



B{r) = B th {r) + B dis (r) 

2(1 - cos(qr)) 



B th (r) = T 
B dis {r) = T 



dq- 



cq* 



Tr 

c 



(52) 
(53) 



(tq 



2(1 - cos(gr)) f 1 dv [a] (v) 



cq z 



o v 2 cq 2 + [a] (v) 

(54) 



The structure factor in -E>dis( r ) is a combination of prop- 
agators ( C(? /[< T K V )+ 1 ) — organized by the RSB parameter 
v, with an increasing self-energy [a] (v) bounded by its 
value at the full- RSB cutoff [a] (v c ). The corresponding 
roughness is eventually computed by integrating explic- 
itly over the Fourier modes q. 



A. GVM roughness of a ID interface 

Using the identity (B16) we obtain an analytical ex- 
pression for a generic [a] (v): 

u n t f 1 dv [in f * 2 ( X ~ cos (<? r )) 
B dis {r) =T • [a] (v) ■ / tiq- 

Ja v 



cq 2 (cq 2 + [a] (v)) 



T ^ (-r/Vc) k ['dv , 



k=2 



(55) 



which simplifies for our full-RSB solution (41) into: 

(r/r ) k ( 1 



k=2 



5k- 6 



+ (1 - Vc) 



(56) 

where ro is a characteristic lengthscale which appears 
naturally in the formalism in order to obtain dimension- 
less quantities. It is defined by: 

[a] (v c ) = cql qa ^° ro = c j [ a ] ( Vc ) (57) 
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and can be made explicit using (41): 

5 5 7T 1 



''0 



T 

3 7 cD 2 [v, 



(58) 



Bdi s ( r ) is thus composed on one hand of the prefactor 
which fixes its dimensions and is actually the thermal 
roughness at the scale ro, and on the other hand of a 
dimensionless series in (r/ro) including the parameter v c 
(45). 

The whole displacement correlation function B(r) is 
plotted in Fig. 3 where we distinguish the low versus 
high temperature cases. In Fig. 4 we show the evolution 
of the roughness at increasing T and fixed D versus the 
other way around. A summary of the different roughness 
regimes along with their corresponding exponent £ and 
their crossover lengthscales is given by Fig. 8(a). 



1. At small lengthscales: thermal regime 

At small lengthscales the linear term of B t ^(r) domi- 
nates the whole roughness, and the ID interface fluctu- 
ates as expected as if there were no disorder: 



B(r) r «° B th (r) = — ~ r 2 ^- 
c 



(59) 



This defines the thermal regime of the roughness, of cor- 
responding thermal exponent fth — 1/2. 



2. At large lengthscales: RM regime 

To make explicit the asymptotic behaviour at large 
lengthscales, we can reformulate the truncated alternated 
series (56) using the following definition of the Euler T- 
function and its generalization: 



I» = / dt- ■ e"' 
Jo 



(60) 



T(a,z) = I dt-t"- 1 ■e"' 



which yield: 



B dis (r) = ^v-i 
c 



(1 - w c ) (e- r/ro - r/ro - l) - r/r 



1 



(r/r ) 



6/5 



6 I 5/ 5 I 5' 1 



(61) 



Since T(— 6/5, a) tends exponentially fast to for increas- 
ing a, at large lengthscales the total roughness is domi- 
nated by the following power-law behaviour: 



B(r) 



B, 



asympt 



_ Tr _! 



6/5 



(62) 



with an 
0.970191, 

4+rn 



overall numerical factor (— |T (— |)j w 
and the Flory roughness exponent Cf = 
,i = 3/5 (see Appendix A). This asymptotic 



exponent which is not the exact one Crm — 2/3, is a 
consequence of the variational approximation. It leads 
in particular to an asymptotic temperature independence 
when rg is made explicit using (58): 



B. 



asympt 



1/5 



. c - 4 / 5 D 2 /5 . r 6/5 



(63) 



We will examine this point in more details in Sec. IV E. 



B. Larkin length, RSB and effective width 

The RS section u > v c corresponds to small length- 
scales (large q) in which the interface fluctuates purely 
thermally, whereas the full-RSB section u < v c corre- 
sponds to large lengthscales (small q) and encodes the 
disorder-induced metastability experienced by the inter- 
face in the RM regime. So the cutoff v c , and consequently 
ro must correspond to the definition of a characteris- 
tic crossover lengthscale between two roughness regimes, 
namely the Larkin length 42 , which marks the beginning 
of the RM regime. 

Using the low- and high-temperature limit of [a] (v c ), 
respectively (49) and (50), analytical expressions of the 
characteristic lengthscale r are accessible at finite c and 



r » 2 



1/3 



£ 5/3 c 2/3 L>~ 1/3 



T->c 



2 5 7T T 5 

~3 2 ~~c~D 2 



(64) 



(65) 



At low temperatures it is temperature-independent and 
depends explicitly on the width £, whereas at high tem- 
perature it is the opposite, and the microscopic parame- 
ter £ is completely washed out. This leads in particular to 
a radical simplification of the prefactor Tl " - ' 
at low temperature: 



^v- 1 inBdisW 



k=2 



1 



5k - 6 



1 



which gives at r = tq: 



B dis (r ) w -e 2 



(66) 



(67) 



It follows from this last expression that ro satisfies the 
usual definition of the Larkin length L c , i.e. the length- 
scale at which the mean displacement of the interface is 
of the order of its width : 



i(L c ) ~ y/B(L c ) ~ £ 



(68) 



At higher temperatures, if we want to keep r Q as the 
Larkin length, i.e. associate the beginning of the RM 
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regime and the breaking of the replica symmetry by im- 
posing ro = L Cl we have then to define an effective width 
£ e g with respect to the roughness at r = r : 



Bdis(r ) ~ 

3 T 3 
2cD 



Tro 
c 



1 



eff 



(69) 
(70) 



so at high temperature the interface 'thickens' due to the 
increasing interplay of thermal fluctuations with disorder. 



C. Crossover at low temperature: modified Larkin 
regime 

The end of the thermal regime can be defined as the 
lengthscale at which the quadratic term ~ r 2 in the sum 
of power-law corrections of Bd± s (r) becomes of the same 
order as the linear term of the thermal roughness: 



2V r . 



n = r 



5/4 - v c 



(71) 



Note that r\ does never diverge since v c < 5/6 < 5/4. 
At low temperature it can be made explicit using (64): 

ri ^(Stt) 1 / 3 ■ T ■ £4/3^/3^-2/3 (?2) 

9 

This increases linearly with the temperature and at some 
point will even merge with rg, defining an upper critical 
temperature with the help of (45): 



n = r Q 



5 

12 



2 9 7T 



1/6 



w 0.42 
(73) 

Up to a factor 1/2 this criterion is equivalent to (48). 

For intermediate lengthscales between r\ and tq the 
roughness function shows a smooth connection between 
the thermal regime at r < r\ and the RM regime at 
r > r"o, as can be seen in the low temperature case 
Fig. 3(a). At fixed strength of disorder D and increasing 
T the crossovers r\ and r Q are getting closer and squeeze 
this intermediate regime, as can be seen in Fig. 4(b). 
In other dimensionalities this intermediate regime would 
be described by the Larkin model 42 , which predicts a 
power-law behaviour of exponent = 4j = 1 ^ . For the ID 
interface (d = m = 1) the perturbative expansion of the 
Larkin model is not valid, and this regime is described 
by a sum of power-law corrections in r, starting from the 



quadratic term 



in -Bdis( r ) (56). By extension we 



thus call it the modified Larkin regime. 

In the zero-temperature limit the linear roughness 
Bth( r ) = c 1S su PP resse d and the small lengthscales 
are dominated by the quadratic term ~ r 2 , which cor- 
responds to an exponent £ = 1 which is again not the 
Larkin exponent Ci|d=l = 3/2. 



D. Crossover at high temperature 

At higher temperatures the modified Larkin regime has 
disappeared, and the intersection of the thermal and the 
asymptotic roughness, respectively (59) and (62), defines 
another characteristic crossover length r*: 



1 



r ■ v c 



■§r(-i/5) 



(74) 



and since the value of the cutoff v c saturates at 5/6 at 
high T, up to a factor 7%/ro s» 0.414225 the crossover 
r* is equivalent to the Larkin length r . Using (65) we 
eventually obtain: 



T-joo 2 T 5 
~ 5 ' ~c~D 2 



(75) 



This result can be predicted by scaling arguments as- 
suming that there is no intermediate regime between the 
thermal and the asymptotic Flory-like regimes. This hy- 
pothesis is true only in the high-temperature regime, e.g. 
imposing £ = from the beginning (see Sec. VI). This 
last approximation actually wrongs the low-temperature 
physics below the Larkin length since it misses the in- 
termediate modified Larkin regime, whereas the GVM 
clearly points out its existence. 



E. Alternative formulation of the GVM roughness 

We finally give an alternative formulation of the GVM 
roughness (56), starting from (51) and (BIO), and rescal- 
ing in the propagators the energy per Fourier mode 
G~ 1 (q) = cq 2 with respect to the mass term at the RSB 
cutoff [a] (v c ). The ad hoc changes of variable make ex- 
plicit the 'crossover function' from the thermal to the RM 
roughness regimes and provide a straightforward way of 
predicting the temperature dependence of the asymptotic 
GVM roughness. 

The following argument being quite generic, we start 
with an internal dimension d and a generalized isotropic 
elasticity indexed by the exponent fi > (/J, = 2 for the 
usual short-range elasticity (4)), which translates into: 

G~ l {q) = c\q\» H ci [u] = lf^q- G~ 1 (q) ■ 

(76) 

As seen in Sec. IIIC, the GVM procedure constrains the 
form of a full-RSB mass term [a] (u) into the following 
power law (y > 0): 



(XL \ v 
7^p) , M (U > Vc) = [(A (v c ) 

(77) 

where v = 10 and x\> = 1 in (37). 

The thermal roughness (53) generalizes into: 

B th (r) = -J 3 d q ~ -q * \ q „ 1/r 

(78) 
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B(r) 



asympt 

(r) 

.6/5 




n~ t 



ro ~ T° r 

RANDOM 

THERMAL MODIFIED ! MANIFOLD 

&. = 1/2 LARKIN ; =3/5 



(a) Low-temperature regime 



Sefl 





-^asympt (?*) 






r 

THERMAL 

Cu, = 1/2 

(b) High-temperat 


* ~ T 5 RANDOM r 
MANIFOLD 

Cf = 3/5 

ure regime 



FIG. 3. GVM prediction for the ID-interface static roughness 
B(r), in log — log representations; the slope of the curves cor- 
responds to 2£(r) as defined by (17) (£ = c = D = 1). (a) At 
low temperature (T = 10 -3 ) an intermediate regime appears 
between the small and large lengthscales regimes, while (b) 
at high temperature (T = 10) no intermediate regime occurs. 
The scalings of the different quantities, including crossover 
lengthscales, are recalled directly on the figures. 



i.e. skipping the angular dependence in the cosine's ar- 
gument: 

7Y2Cth u — d 

B th (r) with U = ^—- (79) 

c 2 

As for the 'disorder' roughness (54), separating the RS 
(u > v c ) and full-RSB (u < v c ) segments of [a] (u) and 
rescaling the energies with respect to [a] (v c ) we can ap- 
ply these two changes of variables: 



eg' 



ro\q\Y 



[a] (v) _ ( v/v, 



i.e. r 



(80) 



[a}(v c ) \\q\»/» 



W v (81) 




T = 1CT 4 
£ = c = l 
D 6 {10 _3 ,10 _2 ,10 _1 ,1} 



(a) At fixed temperature 




15 = 1 

e = c =i 

r e {icr 6 , io~ 5 , icr 4 , ur 3 , ht 2 , ht 1 , 1} 

0.001 o.i 10 T 

(b) At fixed disorder 

FIG. 4. Evolution of the ID-interface static roughness at fi- 
nite width £. (a) At fixed T (low-temperature regime) all the 
roughness curves start from the same _B t h(?~) = — at small 
lengthscales, and behave asymptotically in r 6 ^ 5 with a T- 
independent prefactor increasing with disorder, (b) At fixed 
D all the roughness curves collapse on the same Basympt (r) ~ 
T°r 6 ^ 5 at large lengthscales (see (62)), and start from a mi- 
croscopic thermal roughness which is boosted by increasing 
thermal fluctuations. 



We eventually obtain 
Tv- 1 



Sdi s (r) 



„ 2 Cth 
~'o 

1-v, 



_ 2(1 - cos(qr/r )) 



W + 1 
with the 'crossover function' 

F(x) 



+ \q\-^-F(\q\-^) 



dv- 



1 



(82) 



(83) 



As for the cutoffs in the Fourier-modes integral, they are 
simply redefined with respect to the Larkin length ro and 
are thus harmless at least for the ID-interface case. 

In order to extract the behaviour at large lengthscales 
r, we consider the limit q — > 0, i.e. q — > 0: 

|g|-"/" -> oo =» \q\-^ u ■ F(\q\-^ u ) « (gl"^ (84) 
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Thus, by counting in this limit the scaling of q in (82) we 
can relate v to the asymptotic exponent Casympt predicted 
by GVM: 



q 2^ asymp t 



v 



2(C 



asympt 



Cth) (85) 



We can indeed check that for a ID interface (d = m = 
1) with short-range elasticity (/x = 2) and RB disorder 
(Cf = the GVM Prediction Casympt = Cf = 3/5 

and £th = 1/2 are coherent with v = 10 of (41). 

This link between v and Casympt can actually be used 
to predict the temperature dependence of the asymptotic 
roughness, disregarding the numerical and ^4( Cj _d) contri- 
butions: 

B( r ) r ^°° I r 2<th (r/ro) 2 Caaympt = l r -M/^2 Casympt 

(86) 

Combining (80) and (77) we can make explicit the power- 
law T-dependence of tq and eventually of the roughness: 



B(r) 



„0 . r 2U 



(87) 



We can thus conclude immediately that the asymptotic 
roughness has completely forgotten the interface width 
£ since v c (£) has been washed out, and that ip = 1 as 
imposed by the GVM scheme in (IIIC) leads to a T- 
independent prediction. The elastic constant c has to be 
T-independent to ensure the correct small-lengthscales 
behaviour (78), so the only way to modify this last predic- 
tion would be to choose a T- dependent effective strength 
of disorder D(T). This will be the case in the second 
GVM procedure presented in the next section. 



V. ROUNDED TOY MODEL: GVM, 
ROUGHNESS AND CROSSOVER 
LENGTHSCALES 

The ID interface, that we have tackled up to now as a 
static object, can actually be mapped onto a growing, 
dynamical object, namely the directed polymer (DP), 
which consists in an elastic string starting from a base- 
point y = and getting ahead in 'time' t in a random 
potential V(t,y) 23 . The weight of a trajectory y(t') of 
duration t is given by 



exp • 



(3 / dt' 
Jo 



\(d t ,y(t')) 2 + V(t>,y(t>))] } (88) 



so that through the correspondence (t, y) O (z, x) one 
recovers the weight of an interface of 'length' z = t. The 
roughness at lengthscale z is for instance recovered from 
the statistical properties of the DP at finite time t. See 
also Ref. 43 for another approach. 

Our interest goes to finding a simple approximate de- 
scription of the endpoint y(t) fluctuations, at fixed time 
t. This allows to determine the average of observables 
depending on the sole endpoint position y(t) of the DP 



(such as the roughness ([y(t) — y(0)] 2 )), but of course not 
of observables depending on the whole trajectory. The 
idea is to focus on the effective disorder ruling the ex- 
tremity of a DP of fixed length. 



A. Definition of the rounded toy model 

At fixed disorder V and final time t, we approximate 

the weight Z v (t,y) = J Vy' eT m \ y ' ' v \ of all trajectories 
y'(t) arriving in y at time t through an effective free- 
energy F v (t, y) « \ogZy(t, y) where r] = r)(y) is an 
effective disorder depending only on the arrival position 
y of the DP: 



Fr 1 (t,y) = -y 2 



dy' ■viv') +cte„(t) 



(89) 



The first term ^y 2 is simply the (exact) free-energy in 
the absence of disorder. The i-dependent constant arises 
from normalization along y, while the L y n represents the 
effective potential experienced by the endpoint. Ideally, 
one would infer the distribution of rj from that of V, but 
this task is a priori extremely difficult at finite time. 

However, in the infinite-time limit, it is known 23 that 
the exact correlator ([F(t, y2) — F(t, j/i)] 2 ) is proportional 
to | j/2 — Vi \, as provided e.g. by r\ being a delta-correlated 
white noise. For finite time, numerical evidence 44 sup- 
port that this correlator still behaves as \yi ~ y\\ at 
short distance, and it has been proposed 4546 that tak- 
ing a delta-correlated white noise for r\ in (89) indeed 
provides an approximate but good description of the DP 
endpoint statistics. 

In order to tackle the question of the influence of finite 
disorder correlation length on the properties of an inter- 
face, it is thus very natural to examine in details the cor- 
responding model in which the delta function correlated 
noise r\ is now replaced by a white noise with Gaussian 
correlations of width £ and effective strength D: 



*)(v) = 0, v(y)v{y') = D-R i (y-y')\ D =x (90) 

(see (8) for the definition of the correlator i?|). The limit 

£ —> yields a delta-correlated disorder and corresponds 
to the original so-called 'toy model' 47 . A finite £ actually 
rounds the correlator of the free energy F(t, y) around 

V = o. 

These properties thus define a modified 'rounded' toy 
model, in which the average of an observable O writes 



(o(t)}, 



JdyO(t,y) 



,-0F v (t,y) 



-PF v (t,y) 



(91) 



Although much more simple than the full interface 
model, since the fluctuations are now reduced to those of 
the endpoint instead of the total interface, our rounded 
toy model is still not tractable in an exact way, and we re- 
sort to the GVM approximation in the replica approach. 
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We follow here the route used for a previous modified 
toy model in Ref. 48 (in which was assumed a power-law 
asymptotic decay of the disorder correlator). Having re- 
duced the dimension of the model from the directed poly- 
mer to the toy model, we hope that the GVM method 
catches exact scaling exponents at small and large times, 
together with crossover lengthscales. 

Introducing replicas and averaging over disorder one 
gets the analogue of (12): 



(0(t)) = lim / dy 1 {---)dy n -0{t,y 1 )-e 

n^OJ 

with the effective free energy 



(92) 



RD 

F (t'V) = Vl ~ V min l (Wa. Vb) (93) 



a, 6=1 



where at £ = (i.e. for n a delta-correlated white noise) , 
min; would be the minimum function. At finite £ it ac- 
tually writes: 

miri| (y a , y b ) = - (y a + y b -\y a - Vb ||) (94) 

is the absolute value function, rounded close to the 
origin. In Fourier representation, it reads: 



\y\ i = / dX-e 



-A 2 | 2 2(1 - cos(Ay)) 
A 2 



(95) 



Note that we have not given a definite value to the 
strength D of the effective disorder r){y). To fit the 
infinite-time exact result 23 for ([F(t, y 2 ) — F{t, yi)) 2 ), one 
has to take 



D 



cD 



(96) 



This ensures in particular that the disorder strength D 
of the full model (see equation (2))) has the same dimen- 
sions as the constant D of (96) for the toy model. In 
addition, the expression (96) tells us that the disorder 
contribution to the effective free-energy F„(t,y) (89) de- 
pends on temperature, which is in a way expected since 
it is aimed at representing the contribution of all tra- 
jectories arriving at the endpoint, each trajectory being 
sensitive to thermal fluctuations. 



B. GVM and determination of cr(u) and [a] (u) 

To determine the roughness at all times t we follow the 
lines of Sec. Ill A. The main difference is that we work in 
direct space at fixed t so that there are no Fourier trans- 
formations to consider, and the GVM solution depends 
on time t (i.e. of the lengthscale r in the ID-interface 
formulation) . The equivalent of the trial Hamiltonian is 
provided by the quadratic trial free-energy: 



F (t,y) 



1 

- y 
2 ^ 

a, 6=1 



Va ■ Gj(t) ■ y b 



(97) 



parametrized by a n x n hierarchical matrix whose ele- 
ments are: 



G- b 1 (t) = c/t-5 ab 



o~ ab 



(98) 



The diagonal term c/t corresponds to the elastic part of 
the free-energy, while the off-diagonal elements a ab do 
not depend on t since F is local in time. The varia- 
tional method yields the following disorder-independent 
connected parts 



G~ 1 (t) = c/t 



G c (t) = t/c 



(99) 



thanks to the statistical tilt symmetry as in (28) and 
(29), together with self-consistent saddle point equations 
for the off-diagonal elements 



D 



~e + T[G(t)-G a ^ b {t) 



-1/2 



(100) 



The above is equivalent to the saddle point (32) of the 
full interface but with two important differences: i) there 
is no integration over Fourier modes q since we work in 
a direct space representation; ii) the exponent is —1/2 
instead of —3/2, and thus leads to a different asymptotic 
behaviour. 

As in Sec. IIIB we use the more generic full-RSB for- 
mulation with the mapping parameter u 6 [0,1]: 



a{u) = 



D 



C z +T[G(t)-G{t,u) 



'1/2 



(101) 



Following the same procedure as in Sec. Ill C, we obtain: 



(Ta(u)Y 



(Gc\t) + [a](u)y 2D? 
If we are not on a plateau o~'(u) = 0, we thus obtain 

1/3 



(102) 



Tcr(u) 



2D Z 



(G-^ + M («)) 



2/3 



(103) 



Taking again the derivative d u on this last expression, 
in order to reintroduce a u dependence, and identifying 
cr{u) with respect to (G^ 1 + [a] (it)), we obtain 



a(u) 
[a] (u) 



2 3 D 2 

2 4 - 



3 3 7T 



D 



(u/T) 2 
(u/Tf 



G-\t) 



(104) 
(105) 



where we recognize the same powerlaw structure (u/T) v , 
imposed by the very GVM procedure, as in (37) and (38), 
respectively for [a] (u) and a(u). However the effective 
strength of disorder D given by (96) allows to modify 
this temperature dependence, and thus to circumvent the 
GVM artefact which leads to a T-independent asymp- 
totic roughness for the ID-interface GVM (see Sec. VI). 
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Moreover, contrarily to (38) there is no additive constant 
to a(u), the GVM in direct space provides indeed both 
the solution for o~(u) and [a] (u). In the Fourier repre- 
sentation we obtained first [a] (u) and then had to inte- 
grate it to recover <r(u), which was leading to an addi- 
tive constant. As a consequence the solution a(u) has 
to stick to these power-law functions, with possible steps 
and plateaux inbetween monotonous segments. 

As for the GVM of the ID interface, we look for a full- 
RSB continuous solution, without any additional step. 
In fact, by definition [a] (0) = and consequently this 
imposes the existence of a plateau for u below a first 
cutoff (t) which introduces the lengthscale dependence 
into the solution and later on into the roughness via the 
connected part (99): 



u,(t) = ^Y /3 T./ 



2 \2D 2 J 
[a] (u < u*) = [a] (0) = 



a (u < it*) = cr(0) 



2D 2 c 2 



1/3 



(106) 
(107) 

T- l -r 2lz (108) 



The consistency of (101) for u > imposes the existence 
of a second cutoff u c (£) above which a{u) has a plateau, 
similarly to (43). Indeed, using the inversion formula 
(Bll) for u* < u < u c and taking care of the two plateaux 
in J u dv, the saddle point equation (101) requires again 

that the ^-dependence cancels the u c terms as in (45). 
We thus obtain the following polynomial equation for the 
second full-RSB cutoff it c (£) similarly to (45): 



= A (3/4 - u c ) 



rj-<6 



(109) 



2 4 ' 



2 4 ' 



which is equivalent to substituting (104) and (105) into 
(101) at u = u c . As for the GVM of the ID inter- 
face, it can be shown that the factor 3/4 is closely re- 
lated to the GVM prediction for the asymptotic rough- 
ness exponent, which in the case of the toy model will be 
the exact random manifold exponent for a ID interface: 
u c < 3/4= (2Crm)^ 1 . 

An explicit analytical expression for u c € [0,3/4] can 
be obtained for the two opposite limits for A, which at 
fixed £, D > correspond respectively to T — > and 
T -)• oo: 



A^O 



A— >oo 



3^1/4 T _ 3^1/4 T 3/2 

23/2 ' (1^)1/2 _ 2 3 / 2 ' (£cD)l/2 

3/4 + 0" 



(110) 
(111) 



The polynomial equation for it c (£) (109) leads to the 
analogous definition of a 'critical' temperature with re- 
spect to (48) between the low and high temperature 



regime, using (96): 



A = 1 



_2*_ 

3 3 7T 



1/6 



0.76 



(112) 



and this last constant is of order 1. This criterion shows 
again that the limits of the different parameters, in par- 
ticular T — > and £ — > cannot be exchanged. Moreover 

the scaling between the parameters |r, £>| is exactly 

the same as in (48) for the GVM on the ID interface. 

Using the inversion formulas of hierarchical matrices 
in lim„_i.0j the physically relevant parameters are [a] (u) 
and cr(0), so this full-RSB solution can be summarized 
by: 



1/3 



orsb(O) = 



-2/3 



JRSB 



JRSB 



(u < U*) = [<t] 
(li* < U < U c ) 



RSB ( u >U C ) = 

2 4 - 



(113) 



3 3 7T 



D ■ (u/Ty - c/t 



To be consistent the Ansatz must also satisfy < (t) < 
u c {£) < 3/4, and the collapse of those two cutoffs actually 
defines the lowest time t c for the existence of a full-RSB 
segment in the solution: 



u*(t c ) = u c (i) 



(114) 



This condition is equivalent to [a] (u*(t c )) — [a] (m c (^)) 
and can be reformulated as: 



3 3 tt c 
~2F D 2 



3 4 wcT 3 
~2~6~1> ' 



T 



(115) 



This crossover time can be made explicit in the two limits 
(110) and (111) plugged in (106). This gives respectively: 



t e 2 1 / V/ 4 ■ ^cD- 1 / 2 

A^oo T 3 C 

D 2 



(116) 
(117) 



For smaller lengthscales t < t c the two plateaux merge 
into a RS Ansatz. Following (B5) and (99), we have 
then G(t) — G(t,u) = G c (t) = t/c and the saddle-point 
equation (101) becomes: 



crs(") = ors(0) = ^~[( 2 
Wrs ( u ) = 



c ) (118) 



So the GVM in direct space of the DP toy model yields 
a solution of the saddle-point equation (101) which de- 
pends explicitly on the time t, and consequently on the 
length of the DP; it is RS for t < t c and becomes con- 
tinuously full-RSB with two plateaux for t > t c . This is 
summarized in Fig. 5. The breaking of the replica sym- 
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ctrsbW-c 2 / 3 ^ 2 /^- 1 -*- 2 / 3 



-c/t 



Mrsb ( w ) 



t \, t 



u*(i c ) = it c (f) 



t <t c 



"11 



D f~ 2 TtY 1/2 



(a) 



3/4-- 



u c (A) 



A 



A 



(b) 



FIG. 5. GVM solution for the toy model, (a) Self-energy 
[a] (u) given by (113). At large times (large lengthscales) 
t > t c it is full-RSB with two plateaux which merge at t = t c , 
as indicated by the black arrows, and at small times (small 
lengthscales) t < t c it is RS. (b) Full-RSB cutoff u c as a func- 
tion of A ~ , ^1 „ ~ . obtained by solving the polyno- 

(£-D) 2 (£cD) 2 ' J ° ^ J 

mial equation (109). It starts linearly at A — > and saturates 
to 3/4 at A — s> oo similarly to Fig. 2(b), but D = D(T) will 
modify its low-temperature dependence (110). 



metry at t = t c (£) has the same source as the full-RSB 
cutoff w c (0 in the GVM of the ID interface: it marks the 
transition from a small-lengthscales regime where the DP 
starts by fluctuating thermally and crosses smoothly over 
to the large-lengthscales RM regime where the disorder- 
induced metastability is dominant 48 . In the next sec- 
tion we compute the corresponding roughness, which also 
gives access to the actual crossover from the thermal to 
the RM regime depending on the explicit finite width £. 



C. Computation of the roughness 

The roughness of the ID interface (15) translates 
into the transverse fluctuations of the DP's endpoint 
((y(t) — y(0)) 2 ), and since y(0) = the corresponding 
GVM correlation function is given by: 



W) 2 )=«mG(t) 



(119) 



Using the inversion formula (BIO) we eventually obtain: 
t > t c : lim G(t) = I [ ^] i 4 / 3 - f 2 (120) 



t < t c : lim G(t) = 




-1/2 



(121) 



We can check that this quantity and its derivative are 
continuous, in particular at the junction of the RS and 
full-RSB solutions: 



lim G ns (t c ) = lim G RSB (t c ) = f " " " "' 



ra-vO 



n->0 



3/4 - u c 



(122) 



3T 



d t lim G RS (i)| t=tc = d t lim G RSB (t)| t=tc = — (123) 

n— >0 n— >0 U c 

In order to compare the GVM predictions for the ID 
interface and the rounded toy model, we use (96) to re- 
cover the physical temperature dependence of D and we 
identify: 



r = t c , B{r) = (y(t) 2 )\ t=r 



(124) 



Keeping track of the RS and full-RSB solution we obtain 
the roughness predicted by GVM on the toy model: 



r > r : B RSB (r) = 
r <r : B RS (r) 



3 / 2L> 2 



2 \irc 2 T 2 
Tr D 



1/3 



,4/3 



c cT^/n 



t) 



(125) 



-1/2 



(126) 



with the characteristic lengthscale associated to the full- 
RSB cutoff u c (|) (109): 



ro 



3 3 tt T 5 



(127) 
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This length differs from (58) essentially by the scaling of 
u c . 

We show in Fig. 6 two graphs illustrating the low ver- 
sus high temperature regimes of B{r). A summary of the 
different roughness regimes along with their correspond- 
ing exponent C, and their crossover lengthscales is given 
by Fig. 8(b). 



D. Roughness regimes and crossover lengthscales 

We recover at small lengthscales the same thermal 
regime of exponent Cth = 1/2 as in Sec. IV A 1, whereas 
at large lengthscales the power-law behaviour is directly 
given by (125) which predicts the exact RM exponent 
Crm = 2/3 of a ID interface instead of the Flory ex- 
ponent Cf = 3/5 obtained in Sec. IV A 2. Note that 
Crm = 2/3 coincides with an exact prediction for the 
original toy model (£ = 0) 39 . This leads in particular to 
an asymptotic temperature-dependence ~ T -2 / 3 , which 
makes the B(r) curves intersect with each other upon 
varying T, as illustrated by Fig. 7. 

As in Sec. IV B the lengthscale r = t c corresponds 
to the Larkin length of the system (associated to the 
appearance of full RSB), at least if we define an effective 
width £ c ff at high temperatures. Indeed, the expansions 
(116) and (117) give the following analytical expressions 
for ?*o : 



r « 2 1 /V/ 4 • |3/2 c£r i/2 T i/2 (128) 

(129) 



D 2 



which, when compared to (64) and respectively to (65) 
predicts exactly the same scaling of parameters at high 
temperatures, but introduces a temperature dependence 
~ T 1 / 2 and slightly modifies the scaling-exponents values 
at low temperatures. Combining (122) with (109) the 
roughness at tq is given by: 



~ 2 3/2 + u c 3 3 T 6 
B(r )=£ ■— = — 7T- 



3/4- 



2 4 (cD) 2 



(130) 



so we have respectively, in the two opposite temperature 
limits (110) and (111): 



B(r ) 



T ~° if 2 



, , T^oo 2 4 T 6 ~ 2 
B{ro) * J n (cW 2 =^ 



(131) 
(132) 



Below the Larkin length the roughness is described by 
the RS solution: 



B RS {r) = Bth(r) 



ecD 



^T 3 

r c = c£ 2 IT «■ B th (r c ) = £ 2 



f(r/r c ) 



(133) 



m = 



VTTx 



where the dimensionless f(x) describes the actual shape 
of the crossover along with the definition of r c . The ther- 
mal roughness obtained from the rounded toy model dif- 
fers on two points from its counterpart for the ID inter- 
face, as a consequence of the GVM in direct-space versus 
Fourier representation: i) the RS solution encodes both 
the thermal and the intermediate crossover regimes, in- 
stead of the thermal regime only, ii) the crossover func- 
tion is given explicitly after the retrieval of the thermal 
roughness, whereas in (56) it was given as a series. Taylor 
expanding f(x) around x — 0, we have: 



Blr) fa hr 2 = r 3 =- 



-0(Tr 4 ) (134) 



so the end of the thermal regime can be defined similarly 
to (71) as the quadratic take-off ~ r 2 : 



Tn 2 D 
c ^cT 



T 2 | 



(135) 



The collapse of r\ and tq yields the analogous criterion 
to (73) for a 'critical' temperature: 



n = r a 



37T 1 / 6 T 

2 4 / 3 (|c£>)V3 



(136) 



which is compatible with the previous criterion (112) on 
the GVM solution itself. 

As for the zero-temperature limit, it is ill-defined, since 
r — > and Bdis(?**o) is finite (131), whereas its deriva- 
tive and -BrsbM itself clearly diverge according to (123) 
and (126). The toy model is actually an effective model 
at finite temperature, and is thus non-physical for this 
particular limit. 

Finally at high temperatures there is a single crossover 
lengthscale between the thermal and the RM regime, de- 
noted r* and scaling as (74): 



TV* 



2D 2 

KC 2 T 2 



1/3 



„ 4 /3 



4tt T 5 
27 cD 2 



(137) 



which is equivalent up to a factor 1/27 to the Larkin 
length ?*o at high temperatures and can be predicted by 
scaling arguments. 



Remarkably, except for the numerical factors, both these 
expressions are strictly equivalent to (67) and (69) ob- 
tained in the previous GVM procedure. This is com- 
patible with a generic scaling analysis of the model (see 
Sec. VI). 



VI. SCALING ANALYSIS 

In order to explain the universal expressions of T c , r 
and £ e ffj we follow in spirit the scaling analysis presented 
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B, 
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(a) Low-temperature regime 
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-^asympt ( ?* ) 
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r\ 






r* ~ T 5 
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(b) High-temperature regime 

FIG. 6. GVM prediction for the toy-model static roughness 
B(r), in log — log representations; the slope of the curves cor- 
responds to 2£(r) as defined by (17). It can be compared to 
Fig. 3. (i = c = D = 1, (a) T = 1(T 3 , (b) T = 10) 




FIG. 7. GVM roughness at fixed disorder for the toy model, 
on four orders of magnitude for the temperature: the thick 
curves correspond to -Brsb(V)> the thin ones to Bns(r), and 
they connect at the Larkin length ro. Unlike Fig. 4(b) the 
curves intersect because of the asymptotic scaling ~ y~ 2 /3 
and become thus unphysical in the zero-temperature limit. 



in Refs. 28 and 29, but keeping £ finite. We fully harness 
the scaling relations of the Hamiltonian and the rough- 
ness so as to derive identities between different scaling 
exponents, and also to show how the Flory exponent £p 
arises naturally even if the so-called Flory argument (see 
Appendix A) does not apply. In particular, we account 
for the role of the GVM approximation in those scalings, 
for both the ID interface and the rounded toy model. 

For the sake of generality, we consider a manifold of 
internal dimension d with m transverse components, sub- 
jected to a random-bond potential V. The weight of one 
configuration reads (see Sec. II A): 



Vu exp 



V(u(z),z)]} (138) 



The distribution of V is assumed to be Gaussian with 
correlations 

V(u,z)V(u',z') = D^R{^-)^ d \z~z') (139) 



£ being the microscopic disorder correlation length and 
R a dimensionless function. Let's perform the change of 
variable 

u = au z = bz (140) 

We note that the elastic part of the Hamiltonian rescales 
as: 



ud-2 a 2 



d d z^(V- z u) 2 (141) 



J d d z c -(v z u) 2 = b 

while (in distribution) the disorder part rewrites: 

d d z V{u, z) = b d ' 2 a- m ' 2 I d d z V(u, z) (142) 



where V has microscopic correlation length £/a. Search- 
ing for a unique exponent £ such that a = b 1 * and such 
that the elastic and the disorder parts of the Hamiltonian 
both scale in the same way, we find 



C = Cf = 



4 + m 

Starting from the definition of the roughness 

:P\V\ 



(143) 



B[z) 



VV- 



z v 



and denoting 



(one has £th 



XF = 2(Cf - Cth) 



(144) 



(145) 



2-ci 



in the short-range elastic case, and 
Xf > 0), one eventually obtains that the roughness obeys 
the following scaling relation 

B(z; c, D, T, f) = 6 2Cf B(&- 1 z; c, D, b~ XF T, b^ F () 

(146) 
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Note that this scaling relation also holds in the GVM 
approach for the interface: indeed, the replicated Hamil- 
tonian (13) and thus its quadratic GVM equivalent (18) 
both rescale in the same way as the original Hamiltonian 
upon (140), with a — b^ F as implied by the scaling of the 
elastic contribution G~ 1 {q) = cq 2 . 

In a high-temperature regime the thermal fluctuations 
wash out the presence of £ and thus the roughness is 
expect to be ^-independent. It follows that one can probe 
a scale-invariant behaviour (e.g. at small or large z) of 
the form 

B{z) ~ cte T^z 2 ^ (147) 

where we have isolated the T-dependence of the prefactor 
of z 2< > , defining the thorn exponent b- Choosing b = 
T 1 / XF we obtain from (146) that in each regime where 
the scaling law (147) applies, the following relation holds 

b = (148) 
Xf 

This last relation is always verified in the thermal regime 
(at small lengthscales z), where £ = Cth — which 
from (145) yields b th = ~ as expected. Denoting by (rm 
the roughness exponent at large z, one has in particular 

Wu = (W9) 
Xf 

Note however that any relevance of £- would a priori 
add a T-dependent prefactor to the power-law behav- 
ior (147). In the large z regime, the value of C de- 
pends on the model details for disorder: we note for in- 
stance that (149) is verified in the ID-interface GVM ap- 
proach, where as noted previously the roughness is then 
temperature-independent at large scale (buM = 0) and 
the roughness exponent arising from the computation is 
Crm = (f, whereas the exact value Crm; * = 2/3 predict 
for (149) b = — | which is compatible with the numerical 
simulations of Ref. 29. 

Our previous considerations are based on exact rela- 
tions - forgetting about cutoffs in Fourier modes - arising 
from scaling. Let's now examine the (usually approxi- 
mate) Flory argument and determine how it fails to gives 
the correct (rm for the ID interface. It consists again in 
searching for an exponent C such that a = and such 
that the elastic (141) and the disorder (142) parts of the 
Hamiltonian both scale in the same way, and in assum- 
ing scale invariance i.e. that the roughness is a unique 
power law B(z) ~ a 2 ^ F ~ z 2< ' F at all lengthscales, for 
b = z absorbing all z-dependence through (140) - see also 
Appendix A. We see however that indeed taking b = z 
in (146), one would find Crm = Cf f° r large z if it were 
true that B(l; c, D, z~ XF T, was independent of z 

for large z. This last assumption is however wrong since 
in general b 7^ 0. 

The scaling arguments we have explicited above can 
also be extended as follows: we now search in (146) for a 
and b functions of c, D and T so as to absorb in f3H all 



the dependence in c, D and T. We first notice that the 
random potential scales in distribution as 

V(u,z)=a,- t ?b-%V(u,z) (150) 

where V is a random potential of Gaussian distribution 
with D-independent variance and of microscopic correla- 
tion length I = £/a, or in other words V(u, z)V(u', z') = 
£- m R((u - u')/|) S^(z - z'). Then one checks that 
PH\u,V) = H[u, V]\ c= d=i is fulfilled provided (see also 
Ref. 29) that: 

a — ( c ~ d D d ~ 2 T i ~ d ) ( 4 +™>*f {, = (c~ m _D 2 T 4+m ) < 4 +™>*f 

(151) 

This yields the scaling form 

B{z; c, D, T, = a 2 B 1 (b^z; a" 1 ^ (152) 

where B\ [z\ £) is the roughness at distance z with c = 
D = T = 1 and disorder correlation length £. Let's exam- 
ine the information contained in the scaling form (152), 
which holds for both the interface and the toy model with 
or without GVM, as directly checked. We first note that 
if a -1 ^ is small enough, the existence of a small disorder 
correlation length can be ignored; this defines a charac- 
teristic temperature 

d 2-d xf 

T c = c*=*Dz=*£Tt (153) 

above which the effects of £ should be irrelevant. For the 
ID interface (d — m = 1), one recovers indeed the same 
characteristic temperature T c ~ (S,cD) 3 of (48) and (112) 
respectively for the ID-interface and rounded toy-model 
GVM approaches. 

Besides, in the now well-defined high-temperature 
regime T > T c the roughness should scale as 
B(z;c,D,T) — a 2 Bi(b~ 1 z) independently of £, and dis- 
play two asymptotic regimes with no intermediate regime 
(see Ref. 29 for a numerical study). The Larkin length 
r at which the 'thermal' and 'random manifold' regimes 
connect can be directly inferred from (151-152) and is 
remarkably independent on the actual values of £th and 
Crm- Indeed r is solution of a (6 -1 r ) 2<;th = a(6 _1 r ) 2CRM 
which gives r = b. One recovers the same behaviour 
of (75) for the GVM of the ID interface^nd of (137) 
for GVM of the rounded toy model. Similarly the value 
£ 2 ff = B(r ) of the roughness at this length is indepen- 
dent of both exponents: Coff = a,. This yields £ c ff = ^jj 
as indeed obtained both in the interface (70) and toy 
model (132) GVM approaches, even though they do not 
share the same value for Crm- In the low temperature 
phase T < T c however, having more than one character- 
istic lengthscale, the previous argument does not apply 
and indeed we have observed that the low-temperature 
Larkin length r$ differ between the two GVM predictions 
for the interface and the toy model. 

One has to make an important observation about ap- 
plying to the toy model the scaling arguments exposed at 
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the beginning of this section. The Flory exponent ensur- 
ing that both parts of the Hamiltonian scale identically is 
£f Y = I instead of (143) ((^ y = | as assumed in the for- 
mal argument of Appendix A). Applying blindly (149), 
one would infer that b toy — 0, whereas the thorn expo- 
nent for the toy model has the correct value b toy = — §■ 
To understand this, one has to remark that (146) is true 
if we replace D by D ~ which induces an additional 
T-dependence ensuring b toy = — g- In a nutshell, we just 
pointed that it is the T-dependence of the effective dis- 
order D seen by the DP endpoint in the toy model which 
allows for the Flory argument and the GVM method to 
yield the correct (rm and bptM exponents. Note that nev- 
ertheless (151) and (152) hold without restriction for the 
toy model. 
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VII. DISCUSSION 



(a) ID interface 



A. Low and high temperature regimes 

In the previous sections we have determined under var- 
ied approximation schemes the expression of the rough- 
ness B(r), not only in the asymptotic regimes but also 
on the whole range of scales r. We have kept the disorder 
correlation length £ finite, and it is now time to fit the 
pieces of the puzzle together and analyze its role. 

An important quantity that has come to light is the 
characteristic 'critical' temperature (defined here disre- 
garding numerical factors) 



1/3 



(154) 



for the one-dimensional interface (d = m = 1; see (153) 
for generic d and to). It is the temperature at which 
the intermediate lengthscales tq and n collapse - we re- 
call that ri marks the end of the thermal regime while 
Tq marks the beginning of the RM regime, see Fig. 8. 
Strikingly, the same T c holds for the interface (73) and 
the toy model (136). Moreover, in their corresponding 
variational computation T c separates the two regimes for 
the RSB cutoff in [cr](u), that is, corresponds to having 
the parameter A of order unity, see Eqs. (48) and (112). 
The generic scaling arguments of Sec. VI indicate it is cer- 
tainly no coincidence that these criteria, albeit disparate, 
in fact all lead to the same T c . 

Let's now depict each of the temperature regime, hav- 
ing in mind to compare our two models at hand (the 
results are summarized in Fig. 8). 

(i) The high temperature regime T 3> T c is character- 
ized by the existence of only two roughness regimes, the 
thermal and the RM ones, which connect at a lengthscale 
r * ~ ^752 common to both GVMs (see (75) and (137)). 
We have seen in Sec. VI how this unique ^-independent 
length 7% was imposed by the Hamiltonian generic scal- 
ing properties, independently of the actual values of the 
asymptotic exponent Crm- For r < r*, B(r) ~ Tr is 
thermal, disorder plays no role and thermal fluctuations 
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(b) Toy model 

FIG. 8. Summary of the GVM predictions for the ID-interface 
and the toy-model roughness about their roughness regimes 
and their crossover lengthscales, at low versus high tempera- 
tures. 



enhance roughness. On the other hand, for r > r* the 
roughness scales as B{r) ~ j , 2b r 2CR M) me aning that it 
now decreases as temperature increases (b = — \ by scal- 
ing, see also Refs. 28, 29, and 49) - physically, it means 
that thermal fluctuations inhibit the interface excursions 
to spread too widely in the random potential. We also 
remark that for the toy model, the roughness and thorn 



exponent obtained by GVM are the exact ones (Crm = § , 
b RM = — I ) while for the interface one gets Casympt = § 
and bRM = 0- I n the high temperature regime the 
rounded toy model, coupled to the GVM, thus provides a 
reasonable description of the interface and its crossovers. 
It is a more efficient mapping than to use directly the 
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GVM on the variational Hamiltonian. 

(ii) In the low temperature regime T < T c , £ becomes 
relevant and two lengthscales r and r% now define an in- 
termediate regime separating the thermal and RM ones. 
The GVM has allowed us to compute the actual crossover 
function of the roughness (82) and (83), versus (133). It 
is instructive to compare the value of ro to the Larkin 
length L c obtained in the Larkin model or through FRG. 
Using for instance the notation of Ref. 31 (equation (4.17) 
with the random force strength A(0) = Z?£~ 3 and the 
disorder correlation length rt = £) one has, 



c 2 e 
De- 



1/3 



2 £ £ 



D 



1/3 



(155) 



which matches exactly the low-temperature limit result 
(64) r ~ ^5/3^/3^,-1/3 for the 1D inter f ace . Qn the 

other hand, the toy model exhibits a different scaling 
r ~ ^ 3 / 2 c£)" 1 / 2 r 1 / 2 , see (128). In the low-temperature 
regime, the direct GVM results for the full ID interface 
provide a better picture than the use of the rounded toy 
model. We can thus use whatever mapping is more ap- 
propriate depending on the temperature regime we are 
interested in. 

Let us emphasize that depending on their order, the 
limits £ — > and T — > lead to different physical regimes, 
since T c crucially depends on £. The question of de- 
termining whether experimentally one lies in the low or 
high temperature regime is thus of particular interest. 
Note that on the numerical side, the high-T regime can 
be probed through the directed polymer 29 for which the 
roughness rescales as in (152). 

It is also instructive to inspect the zero temperature 
limit in detail: on the one hand the toy model has 
again the correct Crm = §, but since < the toy 
model also predicts a non-physical divergence of B(r) 
as T — > in the random manifold regime, although 
the roughness should remain finite and not depend on 
T, as predicted e.g. from the FRG zero-temperature 
fixed point 31 . One may conjecture that the behaviour 
B(r) ~ T ,_2 / 3 r 4 ' 3 remains valid for r > r* only for a fi- 
nite range of temperature, below which the scaling would 
become B(r) ~ 7<V 4 / 3 with an effective thorn exponent 
b RM = 0. To describe this behaviour in the toy model 
approach, one would need to introduce a large scale cut- 
off (possibly not being the system size) to prevent the 
roughness to diverge in the limit T — > 0. 



B. Consequences for the dynamics: the quasistatic 
creep regime 

Although satisfactory from a theoretical point of view, 
the determination of the crossover lengths ro and r\ may 
not yield directly observable predictions, since in many 
instances the experimental resolution for B(r) is insuffi- 
cient to probe such small lengthscales. A context however 
where an indirect measurement could be achieved is that 



of interfaces driven by a small force F (compared to the 
'depinning' force F c ). Instead of following a mere linear 
response, the velocity v(F) of the interface is given by the 
'creep law', 22,31,50,51 archetypal of glassy systems 13,52 



v(F) ~ exp 



-pu c 



(156) 



where /U = d 2 2H ^ 2 ^ is the creep exponent, F c — c^/L 2 
is a characteristic depinning force and U c = c£ 2 /L c an 
energy scale (see Ref. 52 for a review). In the previous 
expressions, £ represents the effective width of the in- 
terface, and its value depends on whether one lies in the 
low- or high-temperature regime - which is in general not 
known in practice. A way to discriminate between the 
two regimes is to use U c itself (readily attainable when 
fitting the creep law (156) on numerical or experimental 
data) since its c, D, T-dependence is directly related to 
that of the width £ and of the Larkin length L c . 

In the high temperature regime, taking U c — c£ 2 ff /r* 
directly leads to a linear behaviour U c = T ', a result at 
first sight compatible e.g. with the numerical simulations 
of Ref. 53 where the measured value of U c is indeed pro- 
portional to T. However, the value of F c used in Ref. 53 
is the zero temperature critical depinning force, not scal- 
ing like the characteristic force F c = c£ G ffA 2 = c 2 D 3 T~ 7 . 
The interpretation of the temperature dependence of U c 
observed in numerical models thus remains to be clari- 
fied. 

In the low temperature regime, one has U c = c£ 2 / 'L c 
where £ is the microscopic width of the interface, and L c 
is given by (155). This yields U c — (£c_D)3 = T c which 
is now temperature-independent. In experiments 1,2,54 on 
magnetic DWs, quantities such as the roughness or the 
response to small force have been measured, and a ques- 
tion is to determine whether the low temperature regime 
describes these results. Before examining this, we first 
point out a possible caveat: it is not clear whether the 
interface is at equilibrium or not; though the measured 
roughness (rm ~ 0.66 is compatible with the theoreti- 
cal value Crm = §, it could also be that the interface 
is stuck in a very slowly relaxing state, reached due to 
the applied external field. In that case, the measured de- 
pinning roughness exponent is also compatible with the 
result £ dop R3 0.63 for an interface with harmonic (V z u) 2 
and non-harmonic (V z u) 4 contributions 55 (this value of 
£ dep j g more relevant for experiments than the ill-defined 
result £ dcp sa 1.2 of the purely harmonic elasticity). 

Assuming nevertheless that the observed interface is 
at equilibrium we can examine whether the analysis we 
have put forward applies. A first evidence that experi- 
ments on magnetic domain wall (DW) fall into the low-T 
regime is that the observed value of U/T actually de- 
pends on the combination of c and £ (see e.g. Table I 
in Ref. 2 where U/T is denoted Td ep /T) as opposite to 
high-T behaviour U c = T. In another experiment on 
magnetic DWs 54 , one can estimate T c from the elastic 
constant c ~ 2.4 • 10~ 12 J m _1 , the domain wall width 
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£ ~ 8 nm and the Larkin length L c ~ 40 nm (estimated 
from small force measurements). Using our expressions 
to compute T c , one arrives at T c ~ 325 K which is in- 
deed slightly above the room temperature T ~ 300 K of 
the experiment, suggesting that the interface is indeed in 
the low temperature regime. Further evidence could be 
provided by a direct evaluation of briM which should be 
zero in that case. 



VIII. CONCLUSION 



compass many correct scalings of the interface, except 
for the T — > limit. One could for instance refine the 
toy model by introducing a large-scale cutoff that would 
constrain the roughness to remain finite in this limit. An 
exact rather than GVM solution may also be attainable. 
Finally, given the success of this model in describing the 
static properties, it would be interesting to ascertain how 
much the dynamics of the interface could be fairly ap- 
proximated by the Langevin dynamics corresponding to 
the rounded toy model. 



In this paper we have examined the static roughness of 
a one dimensional interface subjected to a random-bond 
disorder. Using a Gaussian Variational Method we have 
determined the various regimes due to the finite width £ 
of the interface (or equivalently to the finite correlation 
length of the disorder). We obtained that at large tem- 
perature, there is one universal ^-independent crossover 
length from small to large lengthscales regimes, whereas 
at low temperature two ^-dependent crossover lengths de- 
scribe a more complex the intermediate regime. Results 
summarized in Fig. 8(a) and show that the low tempera- 
ture regime is correctly tackled only if one keeps £ finite. 

To discuss the scaling of the different roughness 
regimes, we have compared the results on the original 
model with the ones on a modified 'rounded' toy model, 
which happens to capture the correct random-manifold 
exponent but turns out to be ill-defined in the T — > 
limit. The high temperature results match those of the 
interface, while at low T the scalings of the intermedi- 
ate lengths differ, due to a distinct RM exponent - see 
Fig. 8(b) for a summary. 

We have described the existence of two temperature 
regimes above and below a characteristic temperature 
T c = (£cD)3. On the numerical side simulations 29 
show for the directed polymer clear evidence support- 
ing the existence of a high-temperature regime with 
B{r) ~ t^I^Crm with b = — |. However the observation 
on the experimental side of a negative thorn exponent b 
or of a T-dependent Larkin length remains to be done. 
One first question would be to know whether the temper- 
ature T c is physically relevant or not. We have estimated 
that in creep experiments on magnetic DWs, the value of 
T c is above, but of the order of T, and it would be inter- 
esting clarify the question e.g by a direct measurement 
of b in both regimes. 

The analysis we have presented also raises several in- 
teresting questions and suggests extensions deserving fur- 
ther inquest: one can wonder to what extent the varia- 
tional method induces artefacts in the crossover length- 
scales scalings we put forward, and try to investigate 
this question by use of the FRG, tackling in particular 
the low- and high-temperature regimes - in other words, 
one would need to reconcile the zero-temperature fixed 
point with the observed T-dependence. On the other 
hand, the modified toy model proved to be an attrac- 
tive simpler version of the problem, rich enough to en- 
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Appendix A: 'Flory-like' scaling arguments 

Simple back-of-the-napkin power- counting leads to the 
so-called Flory (or Imry-Ma) estimation for the rough- 
ness exponent: this is what we relate in that appendix, 
for the interface and for its toy model. Assuming that 
at distance z = L the interface presents typical excur- 
sions of extension u = u(L) in the transverse direction, 
we write that the elastic and disorder contributions to 
the Hamiltonian scale according to 33 

U e] [u} = ~ [ dz- {V z u{z)f ~ L d - 2 u 2 (Al) 
H dis [u,V]= [ dz-V {z,u(z)) ~ L d l 2 u- m ' 2 (A2) 



where we have used that V(z,x) V(z', x') = <5 (d) (z - 
z') R^(x — x') ~ L~ d u~ m to determine the scaling of V. 
Assuming that in the RM regime both contributions scale 
in the same way % c \ ~ "Hdis, we obtain 



i(L) ~ L<* , & 



m 



(A3) 



We analyze in Sec. VI why this Flory-like argument does 
not yield the correct RM exponent (e.g. £f = 3/5 instead 
of Crm = 2/3 for the ID interface). It is instructive to 
observe that on the contrary the same argument actually 
works for the toy model of Sec. V. Assuming now that 
the endpoint of the directed polymer y(t) has excursions 
of order Y at time t, one gets the following scaling for 
the contributions to the effective free-energy (89) 



F el (y,t) = ~y 2 ~ t - 1 Y 2 



F dis (y) 



dy'-n^^Y* 



(A4) 
(A5) 
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where we have used that r)(y)r)(y') = D-R?(y—y') ~ Y^ 1 
to determine the scaling of r\. Imposing F e \ ~ F^is one 
obtains 



y ~ ^ 



>toy _ 



(A6) 



Appendix B: Hierarchical matrices 

In this appendix we recall some properties of the hier- 
archical matrices. Both to fix the notations and for con- 
venience for the reader we also give the inversion formu- 
las of replica-symmetric (RS) and full-replica-symmetry- 
breaking (full-RSB) Ansatz directly in the limit n — > 0, 
using extensively the formulas and derivations of Mezard 
and Parisi in Ref. 33 and 35. 

The replica trick constrains the structure of the n x 
n matrix G7, which must be symmetric and equivalent 
under permutation over the replica indices. The matrix 
has to be inverted in the limit n — > 0. A hierarchical 
matrix can be reconstructed by all the permutations of its 
first line, which can be chosen as the reference sequence in 
which the coefficients are classified monotonously. Such 
a generic n x n matrix G and its corresponding inverse 
matrix G _1 (also hierarchical) are thus defined by: 




) 



(Bl) 



G- 1 



G — ( G, G a 



G-l / 1 — 1 



(B2) 



where G aa = G and G m ' = G" 1 Vo. 

This allows the definition of the connected part of these 
matrices, i.e. the sum of the coefficients on any line or 
column G c = ~Y^ a G a b = ^2b^ab, a conserved quantity 
which satisfies for two inverse matrices the relation 



G c ■ G,. — 1 



(B3) 



The simplest case of a hierarchical matrix is the replica- 
symmetric (RS) Ansatz, in which all the off-diagonal co- 
efficients of the matrix are equal: 

CT 1 = (g^G" 1 ) <^> G = (G,G) (B4) 

and in the limit n — > 0: 

c — Lt y i , Lt c 
G- 1 ~ 1 



G 



(Gc 1 )' 



G 



Gc 1 



1 



G-G 
G- 1 



Gc" 1 



(B5) 




(a)fc = 2 1 - 1 = 1 



(b)fc = 2 2 - 1 = 3 



(C)fc: 



63 



FIG. 9. Examples of hierarchical matrices of n x n blocks, 
with increasing integer k of replica-symmetry breaking. Each 
shade of color corresponds to a different value for the coeffi- 
cients of the matrices. 



If the off-diagonal terms count at least two different 
values {<7oj • ■ • jfffcjj we have a replica- symmetry breaking 
(RSB) Ansatz, the integer k counting the number of such 
breakings (see Fig. 9 for generic examples of this struc- 
ture). The integer n being arbitrarily large, in the limit 
k — y oo the monotonous sequence of coefficients on the 
first line of the matrix is more adequately described by a 
monotonous function G(w), depending on a mapping pa- 
rameter u £ [0, 1]. Thus a full RSB hierarchical matrix 
is defined by 

G = (G, G(m)) with u e [0, 1] (B6) 

The full-RSB Ansatz is actually the most generic de- 
scription of a hierarchical matrix (in regards to the limit 
n — > 0), since the replica-symmetric and fc-RSB Ansatz 
can be recovered using step- functions for G(u). The pe- 
culiar symmetries of hierarchical matrices allow to de- 
termine generic inversion formulas directly in the limit 
n — > 33 , once the first line of the matrix is given. There- 
after they have been adapted to the following definition 
of full-RSB hierarchical matrices: 

d- 1 (q)=(G- 1 -a,-a(u))^G(q)=(d(q),G(q,u)) 

(B7) 

o~(u) being defined as a monotonous function on the in- 
terval [0,1], the discrete sums of matrix operations are 
replaced by integrals. Note that the definition of the 
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connected part G c 1 implies that 



a = — / du ■ <j(u) 
h 



(B8) 



We define the following self-energy, illustrated in 
Fig. 10: 

[a] (u) = u- <j{u) - [ dv a(v) (B9) 



This definition implies in particular that [a]' (it) = u 
a'(u). 



a(u c } 




1 u 



FIG. 10. The mass term [a] (u) corresponds to the shaded 
area between a(u) and the a(u) curve, and thus saturates at 
the cutoff u — u c . 



This self-energy acts in fact as a mass term in the prop- 
agators of the following inversion formulas: 



G 



1 



G-G(u) 



(BIO) 

o v 2 G, 1 + [a] (v) Gc 1 ! 

dv 



u Gc L + [a](u) J u v 2 Gc L + [a](v) 

(Bll) 

G-G{u)= 1 + tdv ^ 2 

Gc 1 + [a] (1) J u ( G -i + [a] {v) ) 2 

(B12) 

The relation (B12) can be obtained from (Bll) by a sim- 
ple integration by parts. 

From the full-RSB point of view the RS Ansatz is the 
particular case when a(u) is a constant: 

a{u) = croVw e [0, 1] => [a] (it) = OVu G [0, 1] (B13) 



in which case the previous inversion formulas collapse 
indeed on the RS case (B5). 

In the context of a GVM in a Fourier-space representa- 
tion (as the one described in Sec. Ill) , those inversion for- 
mulas have still to be integrated over the Fourier modes 
q in order to deal with the saddle-point equation or the 
computation of the roughness. Inverting conveniently 
the order of integration over the Fourier modes and over 
the RSB parameter it, and pushing the ultra-violet cutoff 
to oo, the following identities for propagators are useful 
(A > is typically a self-energy [a] (it)) for the case of a 
ID interface: 



dq ■ 



1 



cq 



A 



dq ■ 



A- 1 ' 2 

A -3/2 



(cq 2 + A) 2 4^ 



(B14) 
(B15) 



and 



dq 



2(1 - cos (qr)) 



cq 2 {cq 2 + A) 

_L_ ( e -y/A/c-r _ i + J^T C . r ) 
J Ac \ v ' ) 

s/A/c ■ r\ 



Ay/Ac 

1 

AVAc 



(B16) 



k=2 



at first for J K dq ■ (^G(q) — G(q, u)j in the manipulation 
of the saddle-point equation and secondly for J R dq ■ 
(1 — cos(qr)) ■ G(q) in the computation of the roughness. 

For the sake of completeness here is finally the formula 
for the trace of hierarchical matrices in the limit n — > 0, 
adapted from (AII.ll) of Ref. 33 to our conventions (with 
[G] (it) being defined similarly to [<r] (it) (B9)): 



lim -TrlogG -1 

n->0 n 



log G- 



lim -TrlogG 

n->0 n 



am 

Gc 1 



x du, fGV 1 
— log 



[a] (u) 



Gc 



log(G c ) 



G(0) 
G r 



1 du 



log 



G c - [G] (u) 
Gc 



(B17) 



(B18) 



The quantity ^Tr log G is actually to be used in the GVM 
for the computation of the free energies J-q and T val as- 
sociated to a trial Hamiltonian H , see Appendix C in 
Ref. 32 for an example. 
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